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^ ■ Abstract 
cn 

o 

^f) . We study the net dielectron production rates from an ensemble of thermal 

! mesons, using an effective Lagrangian to model their interaction. The cou- 

I pling between the electromagnetic and the hadronic sectors is done through 

the vector meson dominance approach. For the first time, a complete set 
of light mesons is considered. We include contributions from decays of the 
^ . type V (PS) PS (V) + e"*" e~, where V is a vector meson and PS is a 

pseudoscalar, as well as those from binary reactions PS + PS, V + V, and 
V + PS — > e+e". Direct decays of the type V e+e~ are included and 



X 



' shown to be important. We find that the dielectron invariant mass spectrum 

naturally divides in distinct regions: in the low mass domain the decays from 
vector and pseudoscalar mesons form the dominant contribution. The pion- 
pion annihilation and direct decays then pick up and form the leading signal 
in an invariant mass region that includes the p — uj complex and extends up 
to the (p. Above invariant mass M ~ 1 GeV other two-body reactions take 
over as the prominent mechanisms for lepton pair generation. These facts will 
have quantitative bearing on the eventual identification of the quark-gluon 
plasma. 

PACS numbers: 25.75.+r, 12.38.Mh, 13.75.Lb 
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I. INTRODUCTION 



One of the ultimate goals of high energy heavy ion physics is the formation and ob- 
servation of a quark-gluon plasma (QGP) as predicted by QCD. A vigorous experimental 
program is under way and it is fair to say that this area is one of the most active fields of 
contemporary subatomic physics. The creation of such a novel state of matter represents a 
considerable challenge, both in its experimental realization and also in the theoretical inter- 
pretation of the experimental results. The lifetimes involved are of the order of ~ 10 fm/c 
and the detailed dynamics of the collision process may furthermore play an important role, 
complicating the extraction of a clear signal. Nevertheless, much progress has been made 
both in theory and in experiment, and we may say that even in the absence of a genuine 
QGP the study of hot and dense hadronic systems is still a fascinating subject from which 
a great deal can be learnt. 

For a while now, electromagnetic signals have been known as ideal probes of strongly 
interacting matter at high temperatures and densities This owes to the fact that once 
they are produced, they will travel relatively unscathed from their point of origin to the 
detector. Since production rates are rapidly increasing functions of temperature and density, 
these electromagnetic signals provide valuable information on the hot and dense phases of 
the reaction. It is hoped that, because of these facts, those signals should constitute precious 
aids in the process of analyzing the behaviour of hot quark-gluon matter 0. As with any 
possible experimental signature of the QGP, a great deal of care must go into the calculation 
of a corresponding "purely hadronic" signal, that is a contribution to the same experimental 
observables from sources other than the deconfined, chiral-symmetric phase. As far as the 
quark-gluon plasma is concerned, one may refer to these sources as the "background" . 

In this paper, we are concerned with the thermal rate of dielectron emission only but 
our treatment is completely general. The source is a hot environment of several meson 
species: for the first time, we use a rather complete set of mesons, rather than restricting 
ourselves to the usual pion gas approximation. The equilibrium assumptions inherent to the 
approaches similar to the one been used here have to be carried to their logical conclusion: in 
such scenarios, once the temperature has been set one can clearly calculate the population 
of species present. These mesons can then interact among themselves, or even decay, to 
produce lepton pairs in the final state. It is important to realize that we deliberately make 
no attempt here to connect with experiment because our calculation is rather meant to 
answer a well defined theoretical question: what is the electromagnetic emissivity (in the 
dilepton channel) of a hot hadron gas? To answer this question, we shall proceed along the 
lines of a similar calculation for photon rates |@. 

We estimate the rates of producing lepton pairs using relativistic kinetic theory. The 
mesonic interactions are modelled with an effective Lagrangian and the coupling of radiation 
to hadronic matter is done in the vector meson dominance (VMD) approach. The values of 
the coupling constants involved are adjusted so that the experimentally measured radiative 
decay widths are reproduced. We describe the details of our model in section |T[ We give 
results in section |T| and finally we end with a discussion in section fV[ 
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II. THE MODEL 



Our starting point is an ensemble of mesons in thermal equilibrium. We consider the 
lightest and thus most abundant strange and non-strange mesons together with their main 
interaction channels. This means we shall include: tt, 77, p, cj, 77', 0, if and K* . The charge 
states are not labelled but all of them are present. This collection can be further divided in 
two categories: pseudoscalar (PS) and vector (V) particles. From this hot meson gas, how 
do we calculate what is the lepton pair radiation? It has been shown that the thermal 
production rate for electron-positron pairs is related to the imaginary part of the retarded 
photon self energy by 

1 
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Here p+ and p_ are the positron and electron momenta, = {E, k) is the virtual photon 
momentum, T is the temperature, and we have set the electron mass to zero (nonzero lepton 
mass is easy to include). R is the number of times per unit four- volume an e~^e~ pair of 
invariant mass M is produced with the specified momentum configuration. Note that the 
above equation is perturbative in the electromagnetic interaction only; it is a completely 
non-perturbative expression in the strong interaction. 

Furthermore, we shall make use of the VMD model, which states that the hadronic 
electromagnetic current operator is given by the current-field identity 

= -—m^p^ - — mj0^ - —mlu^. (2.2) 

9p 94> 9u. 

The above expression tells us how the electromagnetic radiation couples to hadronic (in 
our case mesonic) matter: by first coupling to one of the vector mesons with some coupling 
constant. In the above, we have kept the p, uj and fields, but in some cases we shall tacitly 
include also higher vector mesons by using phenomenological form factors inspired by data. 
We further need a model for how the mesons interact among themselves. For this, we shall 
use a simple phenomenological approach, inspired by the chiral properties of low energy 
QCD. Such classes of phenomenological Lagrangians have been quite successful in the past 
in the description of low energy hadronic physics |^. We are explicitly interested in the 
interaction between the different possible combinations of vector {V) and pseudoscalar (yj) 
fields. For reasons that will become clear shortly we restrict our discussion to the following 
interaction Lagrangians 0: 

C}^v^ = 9vv^ e,,^p d'^Vd'^Vf^ , (2.3) 

and 

J^'^l^ = 9v^^ V^ifd^. (2.4) 

In the above, the coupling constants are fitted for each field combination, in a procedure we 
now describe. We have a model for how mesons interact among themselves and how they 
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interact with the electromagnetic field. With this approach, let us study a simple radiative 
process like the decay of a vector meson into a pseudoscalar meson and a photon like e.g. 
uo 7r°7. In this model, the process goes via the cupvr vertex, owing to G parity conservation 
at the strong vertex, and the p° couples to the photon in virtue of the current-field identity. 
This corresponds to the Feynman diagram of Fig. 1. The ratio of coupling constants from 
Eqs. ( |2.3|) and ( p.2|) , in this case gujp-n/gp, is adjusted so that the correct experimental 
radiative decay width r(ti; 7r°7) is obtained. Our Lagrangians are then "calibrated" 
through all the following processes: p — > 717, 

u — > 7r''7, rj'y, rj' p°7, rj' wy, — > 777, — > rj'-y, (p n'^-y. One 

realizes [c.f. Eqs. (|2.2| , |2.3| , |2.4| )] that via this procedure, we can only fix the ratio of strong 
to "electromagnetic" (vector meson-photon) couplings. However it is this very combination 
we shall need for our specific application. 

We now integrate our model for interacting mesons with a dilepton radiation calculation. 
If we keep a calculation of the photon self-energy at the one-loop level an evaluation of its 
imaginary part, as instructed in Eq. ( p.l|) , will yield processes of the type V (PS) PS (V) 
7*, PS + PS ^ 7*, V + PS ^ 7* and V + V ^ 7*. Since such tree-level amphtudes can 
be readily computed and that our general field-theoretic treatment for dilepton emission 
has been shown to agree with relativistic kinetic calculations (up to temperature-dependent 
effects in the form factors, which have been shown to be small 0]) we use the latter ap- 
proach. Finally note that the two-body channels listed above will kinematically dominate 
the contributions of the type V + PS ^ PS + 7*, which we shall neglect. The inclusion of 
such processes would correspond to evaluation of the photon self-energy beyond the one- 
loop level. The first attempt to investigate the role of processes with more than two mesons 
involved has recently been made in [|^. We will return to this point later. 

The basic relativistic kinetic expression for the dilepton production rate from a process 
a + b e~^e~ is well known and can be written down as 

R - \r f "^'P" d^P+ ^'p- f f /O.N 

J 2E„(27r)3 2Efe(27r)3 2E+{2n)^ 2E.{2'Kf ^'^ ^ ' 

X \M\' (27r)4 + , 

Similarly, we may write a rate equation for the decay process a ^ h + e^e^: 

dRa^b + e+e- . r f d^Pa d'^Ph d^p+ d^P- f , fi^ 

/ OF /o„n:^ of /o„A.s of /o„N.^ of /o„^3 J'^ + J''> ^^"^^ 

2\ 



dM^ J 2E„(27r)3 2Ei,i2nf 2E+{2-nf 2E_{2tt] 

x\Mf (27r)^ s^(p^-p^-p^-p_) Si^M^ _ (p^ + 



In the above equations, Af is an overall degeneracy factor dependent upon the specific 
channel and the /'s are Bose-Einstein mean occupation numbers. 

These equations are not suitable for numerical evaluation because of the delta functions. 
However, they can be cast, using standard methods of simplifying phase integrals and the 
spherical symmetry in momentum space, into an appropriate form. The dilepton production 
rate for the process a + 6 — > e"^e~ becomes 



E ^ 



(2vr)^ 



/ dE, PafaiEa) / dE, p,h{E,) (2.7) 

J ma J miy 
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where s = = ml + ml + 2{EaEi, — PaPb^), s+ = {ma + m^Y and s_ = (m^ — m^y. The 
cross section (Jab^e+e-{s)i is obtained by an evaluation of an appropriate Feynman diagram, 
Fig. 2. The muhiple integral in Eq. ( |2.7| ) is evaluated by Monte Carlo methods. While it is 
certainly possible to do some of the integrations analytically, we chose the avenue of keeping 
a relatively transparent integrand and we let the Monte Carlo approach handle the numerical 
complexity. Moreover, this approach allows us to evaluate any desirable differential dilepton 
production rate easily (see, e.g., [Q or |^). 

Similarly, the rate equation for the decay process a ^ h + e^e~ is now: 

dAP = (2^ dM^ L ^'''^ '"^^^""^^ ^'-'^ 

X J^^dx [1 + fbiEb)], 

where Eb = {E^E; + PaPlx)/ma, El = {ml + ml - M2)/(2m„) and dT^ ^ b + e+e-/dM^ is 
the differential decay width into the appropriate channel. In Eqs. ( |2.tj| ) and ( ^.81 ) one notices 
the Bose-Einstein final state enhancement, an in-medium effect. 

We also include the direct decay channels of the form V — >■ e+ e^. As we will show, 
their contributions are non-negligible. This is especially true in the case of p e^e~. One 
can show that for such decays 

dR V _^ p+p- 3 Ft/ _^ p+p- m\ 



rfM2 27r2 N M2 

where 



PCX) 

B{M^) / dE f{E) _ ]vf 2 ^ (2.9) 



B(M^) = p — ^-5^^^ — — r . (2.10) 

^ ' ^ (M2 - miy + (myFtot)2 ^ ^ 

The constant (3 fixes the normalization of the Breit-Wigner probability density function. Its 
value is not important here as it enters also the factor 

N= JdM' (^)' B{M') , (2.11) 

which ensures the correct overall normalization based on the experimental value of the partial 
decay width into the dielectron channel, Fy ^ e+e-- The integral runs over the allowed mass 
range. 

In the above equations, my is the vector meson mass and Ftot is its total decay width. 
For the narrow resonances (a;, 0) the latter is taken constant but the p° width is given its 
proper mass dependence. 



III. RESULTS 

The decay channels considered have already been listed: they are the same radiative 
decay reactions V (PS) PS (V) + 7, as used to fix the couplings constants of our La- 
grangians, with the obvious substitution: 7^7*. The V + PS — > e~^e~ amplitudes can all 
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be obtained from the decay reaction amplitudes by crossing symmetry. We list the entrance 
channels anyway for completeness. They are: uo vr^, p it, (f) ^ 'P Vi ^ v'y v'^ 

p 1]', K*K and K*K. For each of the PS + PS and V + V reaction, we follow the following 
approach: their "bare" amplitude is calculated, squared, and finally multiplied by a form 
factor obtainable from experimental data on e^e~ annihilation. 

The topic of form factors deserves here a short discussion. Of course no information on 
time-like form factors is available through the analysis of meson radiative decays into real 
photons. With respect to this issue, we have followed a simple prescription. The time- 
like electromagnetic form factor of charged pion is experimentally very well known [^] and 
some experimental information exists also about those of both charged and neutral kaons 
| rT| . In our calculations of vr+vr", K~^K', and K^K^ annihilation rates we have used a 
recent parametrization [|I^ of these quantities. The vector mesons annihilation channels 
have been given the same form factors as their corresponding (by strangeness and isospin) 
pseudoscalar counterparts. In the case of decays and V + PS reactions, whenever the G 
parity and isospin conservation laws allowed a coupling only to the p° and its recurrences, 
the charged pion electromagnetic form factor |12[ was used. In the other cases, we have 
stuck with a form factor equivalent to a simple pole corresponding to the lightest permitted 
vector meson. Our way of normalizing coupling constants by means of the radiative decay 
widths leads us to a belief that this conservative choice of form factors does not introduce 
too much uncertainty. We made only one exception from the simple rules sketched above. 
In the case of the reaction p + vr ^ e^e~ the rules would lead to a simple cj-pole. It would 
be a rather bad approximation because the threshold of this reaction lies below the position 
of the 0-resonance, which thus becomes extremely important. We take therefore a two-pole 
formula with the relative weight between the uo and (p contributions same as in the kaon 
isoscalar form factor Fs = {Fk+ + Fk-)/2 jg. 

We have performed our thermal hadronic calculations at three temperatures: 100, 150 
and 200 MeV. We feel that those reflect a range of energies that is somewhat reasonable, by 
current theoretical standards. 

The results for V (PS) — > PS (V) e+e^ at a temperature of 150 MeV are shown on 
Fig. 3. Not all the decays are shown, but only the dominant ones. Coupling constants 
arguments aside, the largest contributions will come from the radiative channels where a 
heavy meson decays into a light one and a lepton pair. This is precisely what is observed on 
Fig. 3. The largest contribution up to invariant mass ~ 0.65 GeV is from — > tt" e~^e~ . 
Over this range, p ^ tt e^e~ represents the next-to-leading contribution and the other 
decays are at least an order of magnitude lower. Above 0.65 GeV invariant mass, the only 
decay with phase space left is 7r° e^e~. Note that the widths for the radiative 

decays of the uj and p^ are comparable and are two orders of magnitude larger than that 
for — *■ 7r*^7 1^. Dalitz decay {e.g. rj —>■ 7e"'"e~) is of higher order in a and can thus 
be neglected. However, this argument alone is not totally convincing as one could imagine 
that the rj could be massively produced at such temperatures. We have therefore performed 
a calculation of the contribution from eta Dalitz decay to thermal electron pair yield, using 
the VMD prescription for dT/dM^ fl^ with updated coupling constants. We have found it 
in fact to be orders of magnitude smaller than the channels discussed above. 

For the pseudoscalar-pseudoscalar reactions, on Fig. 4 we display a plot of all the 
contributions, again at T = 150 MeV. The different contributions add up to a signal in 
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which the only apparent structures are associated with the p(770) and the 0, with a shght 
shoulder at the p(2150). The peak in the pion form factor at the p(1700) is washed out by 
the kaon contributions. 

We show the V + V contributions on Fig. 5. Above threshold, the sum of these processes 
outshine the PS + PS ones by roughly an order of magnitude. The structure at M = 2.15 
GeV owes to the corresponding excitation of the p. 

The V + PS reactions are quite numerous, we show the brighter dilepton sources on Fig. 
6, again for T = 150 MeV. The dominant channels are uj + vr'', p + vr and p° + rj. The kaon 
channels are not shown but are roughly the size of the vr + p contribution. The strongest 
signal is from u; + vr*', over the entire invariant mass range considered here. Recall from our 
discussion of the decays that the radiative decay widths of the p and the uj are quite large. 

Finally, the total rate corresponding to the sum of all processes discussed so far is shown 
on Fig. 1^, along with a curve representing the tt+tt" contribution only. We also show the 
net direct decay contribution, summing p e~*"e~, uo e^e~ and e^e~ . The 
radiation from these channels turns out to be quite important. The signal from the decay 
reaction p e^e~ closely resembles the pion annihilation spectrum, which in retrospect 
is quite reasonable. 

In all cases (decays, PS + PS, V + V, V + PS) our findings at T = 100 and 200 MeV 
are qualitatively similar, with a global shift in the rate. For these temperatures we therefore 
present only the total rates (see Fig. H). 



IV. DISCUSSION 

Up to now, thermal calculations of the variety discussed in this paper have rarely gone 
beyond a pure pion gas approximation, usually concentrating on the annihilation channel 
The contribution from thermal meson decays has been considered previously [jl5 



To our knowledge it is the first time that extensive mesonic reactions have been included, 
together with direct decays. 

Comparing the individual contributions from different processes (Figs. 3-6) to the total 
dilepton rate (Fig. |^ one sees that the dilepton invariant mass spectrum naturally divides 
in several parts. At low masses, the decay channels clearly dominate the entire spectrum. 
The crossover to the pion-pion annihilation and direct decay signal occurs just above 0.5 
GeV (at the lower temperature, T = 100 MeV, this crossing point is shifted closer to the 
two-pion threshold). Already at M ^ 1 GeV, the total rate dominates over the pion gas 
approximation result by an approximate factor of 3. At M = 1.5 GeV, those rates differ 
by a little more than an order of magnitude. The difference increases with larger invariant 
masses. One also sees that the net rate at the vector meson positions is also larger than 
in the straight 7r+ — tt^ scenario, owing principally to direct decays and also form factor 
effects. Probably the most striking conclusion of our work is that the "usual" pion results 
for lepton pair production calculation holds rather poorly over all regions of invariant masses 
considered in this work. This statement is true for all temperatures studied here. 

Thus, the rate for M ^ 1 GeV is approximately one order of magnitude larger in our 
calculation than in "conventional" meson background calculations. This enhancement is also 
present in the momentum structure of the lepton signal: Fig. ^ is a plot of Ed^R/d^p for 
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lepton pair invariant masses between 1.1 and 3 GeV. These findings should have important 
imphcations in connection with the plasma signal identification. The conventional window 
for thermal lepton pairs of plasma origin is < M < mj/^ 0, precisely the range 
discussed here. An observation of a signal from an exotic source can only be claimed if 
all other sources are under control. Here these would be identified with the Drell-Yan 
mechanism, open charm decay |jT6| and the thermal background we have considered in this 
work. However, before any more quantitative statements can be made, it is imperative to 
complement our calculations with a dynamical model of some sort, in order to make contact 
with genuine observables. Work in this direction is in progress. 

It is of interest to compare the rates obtained with other similar calculations. Some 
recent interest has been devoted to the emission of lepton pairs from pionic bremsstrahlung 
processes |T^ . It was concluded that the radiation from the external pion lines in pion-pion 



collisions would be a dominant contribution to the low mass lepton spectrum. Comparing 
with pion bremsstrahlung calculations at T = 150 MeV, we realize that the low mass sig- 
nal is the same magnitude as the net meson decay contribution. Correcting the pion-pion 
bremsstrahlung rate for the Landau-Pomeranchuk effect |T^ will cut this pion signal by some 
factor. This factor is only f» 2 for low invariant masses and T = 150 MeV ||19|. This correc- 
tion also goes down as invariant mass grows. This will then leave the pion bremsstrahlung 
to compete with the decay channels contribution, up to the two-pion annihilation threshold. 

In this inquiry, we have pursued the same goals as a similar photon production calculation 
0. Our main aim has been to identify the most important dilepton production processes 
which operate in a hadron gas. We considered only the decays and reactions with the 
minimal possible number of hadrons: one in the decay final states, none in the final states 
of two-initial-hadron reactions. These processes are believed, on the basis of the phase- 
space and order-of-interaction arguments, to be dominant here. The reactions of this kind 
(2^0 hadrons) do not operate in real photon production due to restrictions from energy- 
momentum conservation. However, the dominant reactions for photon production a + 6 — *■ 
c + J can produce virtual photons as well. It is clear that they would populate preferably 
the low-mass region. Even there they would be probably negligible, as it was shown for the 
case of TT + TT 71+ dilepton in |0. But one cannot exclude surprises. The latter process 
amplifies, together with the three pion annihilation channel, the omega peak in dilepton 
spectrum. This may in turn serve as a signature of a hadron gas creation [0]. It has also 
been pointed out that the Ai meson could have a significant infiuence on the real photon 



yield |2^, through the process tt p Ai vr 7. This conjecture has been carefully 
analyzed in a recent paper [^. The refiection in the thermal dilepton sector is certainly 
worth studying as well. Three body initial state processes a + b + c e^e~ may also 
contribute significantly in the high invariant mass region |^. We intend to study all these 
points in detail in future work. 
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FIGURES 



FIG. 1. The Feynman diagram for the radiative decay u tt 7 in the model described in the 
text. 

FIG. 2. A "generic" two-body ampHtudc with lepton pairs in the final state. The different 
{a, 5} combinations we consider are enumerated in the text. The vector meson V is chosen through 
isospin and G parity arguments. 

FIG. 3. Differential rate for lepton pair production via vector or pseudoscalar meson decay. 
The dashed line represents the contribution from u — Tr'^e+e", the dashed-dotted line is the 

rate from p vre+e". The dotted line is the process (p 7r^e^e~. The structure in the latter 
channel is due to the p{770). The solid line is the sum of all the decay processes, including those 
not listed in this caption but enumerated in the main text. 

FIG. 4. Rate from PS + PS type reactions. The dashed line is the rate for the pion annihilation 
process. The dotted curve represents the contribution from K'^ + K~ . The dashed-dotted line 
is the rate from K'^ annihilation. The solid line is the sum of the PS + PS processes. 

FIG. 5. Rate from V + V type reactions. The dashed curve is the p+ + p~ contribution. The 

dashed-dotted and dotted curves represent charged and neutral K* annihilation, respectively. The 
solid curve is the sum of the V + V contributions. 

FIG. 6. Rate from V + PS type reactions. The dashed curve is the rate from co + tt^. The 
dashed-dotted curve is the contribution from p + tt and the dotted curve is the rate from p^ + rj. 
Again, the solid line is a sum of all V + PS processes, as enumerated in the text. 

FIG. 7. The solid line is the total rate at T = 150 MeV from all processes discussed in the 
text. The dashed line is the pion-pion annihilation contribution only. The short-dashed curve 
represents the contribution from direct vector meson decays. 

FIG. 8. Same caption as Fig. 7 but for the temperatures T = 100 MeV (lower curves) and 
T = 200 MeV (higher curves). 

FIG. 9. The lepton pair momentum spectrum, Ed^R/d^p for lepton pair invariant masses 
between 1.1 and 3 GeV. This lower bound is chosen so as to exclude the (f) peak. The full curve 
represents the contribution from all processes described in this work. The dashed curve is the 
pion-pion annihilation contribution only. 
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